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Abstract 

The electron-screening acceleration of laboratory fusion reactions at astro- 
physical energies is an unsolved problem of great importance to astrophysics. 
That effect is modeled here by considering the fusion of hydrogen-like atoms 
whose electron probability density is used in Poisson's equation in order to de- 
rive the corresponding screened Coulomb potential energy. That way atomic 
excitations and deformations of the fusing atoms can be taken into account. 
Those potentials are then treated semiclassically in order to obtain the screen- 
ing (accelerating) factor of the reaction. By means of the proposed model the 
effect of a superstrong magnetic field on laboratory Hydrogen fusion reac- 
tions is investigated here for the first time showing that, despite the consid- 
erable increase in the cross section of the dd reaction, the pp reaction is still 
too slow to justify experimentation. The proposed model is finally applied 
on the H 2 (d,p) H s fusion reaction describing satisfactorily the experimental 
data although some ambiguity remains regarding the molecular nature of the 
deuteron target. Notably, the present method gives a sufficiently high screen- 
ing energy for Hydrogen fusion reactions so that the take-away energy of the 
spectator nucleus can also be taken into account. 
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I. INTRODUCTION 



At astrophysical energies of a few keV corresponding to stellar temperatures of several 
millions degrees kelvin the cross section a (E) of the predominant s-wave fusion reactions is 
given by 

a(E) = ^P(E) (1) 

where the astrophysical factor S (E) embodies all the nuclear effects of the reaction and for 
non-resonant cases is a slowly varying function of the center-of-mass energy E. On the other 
hand , the penetrability factor P (E) embodies all atomic effects of the reaction and when 
the electron cloud around the fusing nuclei is ignored it is given by P (E) = exp (— 2im) 
where n is the Sommerfeld parameter. 

As the astrophysical factor varies slowly with energy we usually replace it with a trun- 
cated Taylor series which will be studied extensively in the present paper 

S(E) = S (0) + S' {0) E + 0.5S" (0) E 2 (2) 

Any error in the zero-energy astrophysical factor S (0) is actually an error in the corre- 
sponding reaction rate in the stellar plasma, which in turn reflects linearly on the energy 
production rate. 

In the past years there have been exhaustive efforts to extend measurements of the S (E) 
towards even lower energies PJ @ in order to obtain a reliable value for S (0). This is neces- 
sary as extrapolating higher energy data to zero energies introduces an inevitable numerical 
error. However, at such low energies, the electron cloud that screens the fusing nuclei en- 
hances the fusion reaction by lowering the Coulomb barrier. Consequently, disregarding its 
presence leads to an overestimation of 5(0). Unfortunately, even very recent experiments 
|| cannot explain the screening enhancement which exceeds all the available theoretical 
predictions as was recently admitted 0,0. 

Various authors have studied the influence of the atomic cloud on the cross section of low 
energy nuclear reaction. A qualitative study 0, which parametrized various atomic pro- 
cesses such as molecular formation, excitation and ionization, yielded a fair approximation 
for the possible contributions of the electronic degrees of freedom in the nuclear collision 
experiment. Moreover, by assuming a constant charge density around the target nucleus, 
a subsequent model predicted a screening shift which was compatible with the experi- 
mental data . However that assumption is an oversimplification which will be amended in 
the present paper. The most sophisticated approach has been a few-body treatment |8j 
which established a lower (sudden) and a higher (adiabatic) limit for the screening energy 



transferred into the relative nuclear motion. Although more studies followed |9j [ID], which 
also extended the calculations to molecular fusion reactions Jll] , despite their mathematical 
rigor they could not explain the discrepancy between experimental and theoretical screening 
energies. 

In this work there is presented a mean-field model for the study of screened nuclear 
reactions at astrophysical energies in the laboratory [|IIJ. That model agrees well with the 
available experimental data, thus enabling us to improve the accuracy of the associated 



2 



astrophysical factor. Moreover, by means of the proposed model the effect of a superstrong 
magnetic field on laboratory Hydrogen fusion reactions is also investigated for the first time, 
yielding the associated magnetic accelerating factor. Notably, the present method gives a 
sufficiently high screening energy for Hydrogen fusion reactions so that the spectator nucleus 
take- away energy can also be taken into account. 



II. SCREENED COULOMB POTENTIALS 

After the pioneering work || that established the importance of atomic effects in low 
energy nuclear reactions various authors have tried to create models that account for the 
observed enhancement. A simple model J7|, suggested at an early stage, assumed that 
the electronic charge density around the target nucleus is constant, thus predicting for the 
nucleus-atom reaction between the atomic target Zieand the projectile Z 2 e a screening 
energy U e = (3/2) ZiZ 2 e 2 a~ 1 . In order to take into account the dependence of the screening 
radius on the charge state of the reaction participants, that model used a screening radius 



0.8853a (z 2/3 + Z 2/3 ) 1/2 (3) 



taken from scattering experiments [T3| so that 

3 , v 2/: 

where the Bohr radius. Although that screening energy is larger than the one predicted 
by the simple formula || U e — ZiZ 2 e 2 (a /Zi) _1 it has some very obvious defects. The 
assumption that the charge density is constant leads to an unnaturally sharp cut-off at a 
distance r = a from the center of the target nuclei, which is not born out either by theory or 
experiment. Moreover, atomic excitations and deformations of the target atom are totally 
disregarded. On the other hand normalizing the charge distribution so that the total charge 
is — Zie gives a charge density 

Po = "J— 3 ( 4 ) 

4 7rcr 

In order to assess the validity of that density we can consider the hydrogen-like atom Z\e 
which will also be used in this section . The charge density at the center of the cloud of such 
an atom (when the electron is in its ground state) is p$ = — e (Zi/a ) 3 /it. It is obvious that 
for Zi = Z 2 = 1 we obtain p ~ — (e/al) and p$ = — (e/ap) /vr, that is the simplified model 
in question overestimates the central density by a factor of 7r. 

Consequently it is obvious that if low energy nuclear reactions are to be treated by means 
of a mean-field potential a more sophisticated treatment is necessary. 

As a first step we consider a more plausible charge distribution: 

P (r) = Po (l - ^) (5) 

which takes into account the depletion of charge with respect to distance from the center. 
The radius a is the screening radius given by Eq. @ and the charge density p at the center 
of the cloud can be found by means of the normalization condition : 
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pa 

/ p (r) 47ir 4 dr = —Z\e 
Jo 



This integral yields a central value of 



Po 



15Z ie 



7TGT 



(6) 



(7) 



Note that for a collision Z\ = Z2 = 1 we have a central charge density po = 7.68 (e/ag) /tt 
which gives an even larger core density than the constant density assumption. An alternative 
approach would be to consider the value p equal to the corresponding hydrogen-like one 
and then calculate the screening radius using Eq. (||) • The latter treatment gives a screening 
radius 
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1/3 



8Zfn 



a 



gives a 



which is independent of the charge of the projectile. For hydrogen isotopes Eq. 
radius of a = 0.842a 

We can calculate the electrostatic energy by solving the equation of Poisson for the above 
charge distribution with the appropriate boundary conditions, so that 



$ ( r ) 



IhZxe 
Y2~ 



3 
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(9) 



Whenever a bare nucleus Z 2 e impinges on the target nuclei surrounded by the electron cloud 
of Eq. (H]) the total interaction potential in the atom-nucleus reaction channel is 



V(r) 



Z 1 Z 2 e 2 15 Z x Z 2 e 2 
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r\ 2 3 fr 
2 \a) + 10 U 



(10) 



Although the above potential energy is more plausible than the constant charge density one, 
a more reliable charge distribution should be considered which could account for various 
other atomic effects as well as for the atom-atom reaction channel. 

Let us consider a hydrogen-like atom with atomic number Z\. When the wave function 
of the electron is given by \l/ n ; (r, 6) then the charge density around the point-like nucleus is 



-e|#„j (r, 



by which it is obvious that both the previous screening model and that of Ref. are 
imperfect. If we solve the equation of Poisson for hydrogen atoms (or hydrogen-like ions) 
whose electron is in its ground (Is) state we obtain 



$00 (r) 

where the screening radius is 



C 6 / T 

— + - 1 + — ) exp (-r/r ) 
r r V 2r / 



(12) 
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(13) 

If a positive projectile Z 2 e interacts with the above screened nucleus then the total potential 
energy is 

Vbo H = + 1 + 7T- exp 14 

r r r \ 2r / V r / 

On the other hand if we assume that the electron is in an excited state (2s) then the 
potential energy is found to be: 

Z x Z 2 e 2 Z 2 e 2 Z 2 e 2 /, 3 r r 2 r 3 \ ( r 



V 10 (r) = i_ + _J_ i + + + exp(-— (15) 

r r r \ 8r o l" r o 64rg / V 2r / 

It should be emphasized that in the derivation of the above potentials we have assumed 
an unperturbed wavefunction of the target nuclei, throughout the tunnelling process. In 
fact at astrophysical energies the electron cloud responds rapidly and by the time tunneling 
begins the nuclei are so close that the wavefunction is actually that of a hydrogen-like atom 
with charge Z\ = (Z% + Z 2 ) and a screening radius Tq = ao/2Z{ 



III. NUCLEAR REACTIONS AT ASTROPHYSICAL ENERGIES 



At astrophysical energies reactions between light nuclei take place via s-interactions, thus 
enabling us to investigate them by means of the WKB. 

If we assume that a bare nucleus Z 2 e collides at very low energy E with a screened 
nucleus whose electron is in its ground state then the tunneling probability according to the 
WKB method is: 



P (E) = exp 



R 



V 00 (r) - Edr 



(16) 



We can assume that the lower limit of the WKB integral is given in terms of the mass 
number A of the reacting nuclei : R = 1.4 (^A^ 3 + A 2 ^ 3 ^j . For most practical purposes this 
lower bound is set equal to zero as all the nuclear effects of the fusion reaction are included 
in the cross section factor. 

The classical turning point can be obtained by equating the relative collision energy E 
with the potential energy of the interaction. The collision energy is set equal to the Gamow 
peak of the corresponding reaction in the plasma so that: 



Vbo (r c ) = 1-220 ■ [Z 2 Z 2 AT 2 ) 1/3 keV 



(17) 



where A the reduced mass number and T 6 the temperature in million degrees kelvin. For a 
wide range of light nuclei we have performed extensive numerical solutions for Eq. flI7|) as 
well as numerical integrations of Eq. (|T6| ) . At astrophysical energies, just as is the case with 
the Debye-Hiickel model in plasma conditions |14]], the results indicate that throughout the 
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potential barrier the potential energy V 00 (r) of Eq. ( [H| ) can be safely replaced by the much 
simpler formula: 



Z x Z 2 e 2 Z\Z 2 e 2 



Voo (r) 



'00 

r a 

Therefore the WKB penetration factor can be written as: 



P (E) = exp 



2y/2j2 rc(E) \Z x Z 2 e 2 Z{Z 2 e 2 



h Jr 



Edr 



a 



(19) 



The equation for the classical turning point is modified accordingly: 



Z x Z 2 e 2 



1.220- [Z 2 Z 2 AT 2 y keV 



1/3 



(20) 



where we have ignored the screening shift given by: 

Z*Z 2 e 2 



«o 



(21) 



It is now obvious that the relative energy of the reaction has been increased by U e . In that 



case the penetration factor can be easily found to be |10 

f ls (E) ~ exp 



™ (E) § 



(22) 



where the subscripts indicate the excitation state of the target atom. If we follow the same 
methodology for the 2s state we obtain 



irn (E) — 
y AE 



f 2s (E) ~ exp 

The much simpler potential model of Eq. ( |10D gives a screening factor: 

fo (E) ~ exp 

with an energy shift of 



(23) 



™ (E) § 



(24) 



\hZ x Z 2 e 2 



(25) 



where a is given either from Eq. (Rf) or Eq. 
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IV. MAGNETICALLY CATALYZED SCREENING 



By now it is obvious that any shift U e « E of the interaction potential energy V (r) 

V(r) = ^^-U e (26) 
r 

accelerates the fusion cross section of hydrogen isotopes by a factor / ls (E) given by Eq. 
(|22D . That observation will prove very useful in the study of the effects of a superstrong 
magnetic field on laboratory hydrogen fusion reactions which follows. 

As a matter of fact under such extreme conditions the electron-screening cloud is de- 
formed in the sense that it becomes compressed perpendicular and parallel to the magnetic 
field so that the screening potential energy for the strongly magnetized hydrogen atom is 
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e 2 x poo exp 



2 \ l+u ~ a 2 +u J 



p V27T Jo (1 + u) V« + u 

where p, z are the coordinates in a cylindrical frame of reference whose origin coincides with 
the point-like nucleus of the hydrogen atom. 

The natural length unit in the above formula is of course the cyclotron radius so that 
p = p/p,z = z/p , and a is a parameter which depends on the magnetic field and is 
determined by the variational method. The above formula was shown to be reliable for very 
strong fields whereas it becomes inaccurate below the threshold of the intense magnetic field 
regime given by: 

B IMF = 4.7 x 10 9 G. (28) 
In Ref. [15] potential ( |2"7D was applied at zero relative energies in order to obtain the mean- 



life times of hydrogen isotopes in neutron star surfaces. However, a more recent work JT6 
used that potential in a problem where the relative energies were of the order of keV showing 
that for energies E > 0.5keV and fields of the order of B 12 = 0.047 (Bi 2 being the field 
measured in 10 12 G) the classical turning point is so deep inside the cloud that the screening 
shift can be considered constant and equal to the value of the potential at the center of the 



cloud given in Ref. [15 



e 2 2 hi f a + Va 2 - l) 
U e (0,0;a) = ^ = V r^-^ - (29) 



P V27T V« 

In the present work that approximation has been tested for various other fields and energies. 
The results show that for fields as high as B\ 2 = 4.7 and interaction energies E > 0.5 keV 
the screening effect is independent of the angle at which the projectile enters the electron 
cloud and can be considered equal to Eq. ( p9[) . 

Therefore if the target hydrogen nuclei are in such a magnetic field the reaction is going 
to be accelerated by a factor 
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f u (E) ~ exp 



nn (E) 



u e (Q, Q; a) 

E 



(30) 



Figures 1 and 2 depict the acceleration of the pp and dd reactions respectively for various 
magnetic fields and interaction energies. Especially for the pp reaction it is obvious that 
even in such a strong field the cross section is still significantly small. Namely, as the corre- 
sponding zero energy astrophysical factor is S PP (0) ~ 4 x 10~ 22 keV — barns, the screening 
effect in a superstrong field B\ 2 = 4.7 can only increase S pp (0) by roughly one order of 
magnitude compared to the unmagnetized case. 

The dd reaction, on the other hand, can be significantly affected by such a magnetic field 
as it is already much faster than the pp one. At very low energies the increase can be as 
high as two orders of magnitude compared to the unmagnetized case. 



V. THE ASTROPHYSICAL FACTOR OF D-D NUCLEAR REACTIONS. 

Despite the fact that the reactions H 2 (d,p) H 3 , H 2 (d,n) He 3 have been investigated 
since the early days of accelerators [17]] |l8j |pJj , the effect of screening on the associated 
astrophysical S (E) , which will eventually be used in theoretical calculations, is still under 
investigation. In the discussion that follows we will show that our model is compatible with 
the experimental data of that reaction. 

The appropriate treatment of a low-energy experiment should take into account screening 
effects in order to calculate the respective values of S (E) . As a matter of fact once a 
screening model and the associated screening energy U e are adopted the corrected bare- 
nucleus astrophysical factor of the experiment is actually given by 

S b (E) = Eo (E) exp (27m) exp (31) 



Then Eq. (||D is fitted to the data corrected through Eq. (31) in order to obtain the zero- 
energy coefficient S (0) . 

Any effort to extrapolate from higher-energy data or fit all the uncorrected data with 
formula @ is bound to induce errors. 

There are three different ways to analyze low energy fusion data which must of course 
be consistent with each other. We will apply those methods on the available data |2(| for dd 



reactions (E > 2keV) and compare them with the analytic model proposed in the present 
paper. First we note that for energies E > 20keV any screening correction is meaningless 
since the exponential term of Eq. ([31]) is very close to unity at such high energies. Therefore 
we can obtain the asymptotic behavior of the astrophysical factor by using the available 
high-precision experimental data pTI] for higher energies which yielded 



S b (E) = 55.49 (0.46) + 0.094 (0.0054) E (32) 

We can now reasonably assume that this should be a fair approximation of the bare-nucleus 
astrophysical provided its use consistently describes the low-energy experimental data. In 
fact the screened value of S (E) will now be given by 
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/ U as \ 

S (E) = (55.49 + 0.094E) exp (j™-?-) (33) 



where the screening energy U£ s is determined by fitting Eq. ( |33"D to the uncorrected data of 
Ref. |3, so that V a e s = 0.019 (0.003) keV with X 2 = 0.028. 



The second method which will corroborate the validity of the proposed models entails 
fitting all four parameters S (0) , S (0) ,S (0) , U e simultaneously to the uncorrected exper- 
imental data. Thus we obtain a screening energy of Iff 1 = 0.017 (0.003) keVemd a bare 
nucleus astrophysical factor: 

S b (E) = 54.54 (1.39) + 0.608 (0.265) E - 0.026 (0.026) (34) 

with x 2 — 0.011. Obviously, the two previous approaches give results which are compatible 
with each other as expected. Figure 3 shows that both the previous two fits provide a 
satisfactory description of the screening effect. 

The third method is a straightforward application of the theoretical models derived in 
the present paper. However, in order to apply those models on the experimental data 
we have to take into account that the data refer to a molecular target while our models 
refer to atomic ones. Hence, we have to allow for the energy which will be carried away 
by the spectator nuclei plus the reduction due to the molecular binding energy. Although 
this assumption has been argued against |i~lfl , the actual energy reduction for a deuteron 
molecular target has been calculated |^3| by a Coulomb explosion process to be of the order of 



44 eV. Therefore modifying our models for a molecular deuteron target we derive a screening 
energy U e = O.OlO/ceV^ (Eq. fl2"ip) and U e = 0.016 keV (Eq. (^51)) which are in reasonably 
good agreement with the experimentally obtained values. We can now fit the formula 

S (E) = [S (0) + S' (0) E + 0.5S" (0) E 2 ] exp Un^j (35) 

by using the screening shift of our models. The results are as follows 
U e = 0.010 

S b (E) = 57.3 (0.41) + 0.160 (0.125) E - 0.0056 (0.002) E 2 (36) 

withjc 2 = 0.013 and 
U e = 0.016 

S b (E) = 54.93 (0.38) + 0.537 (0.1149) E - 0.0225 (0.007) E 2 (37) 
with x 2 = 0.011 

Although our models are fairly compatible with the experiment there is an inevitably 
degree of uncertainty in the associated astrophysical factors due to the actual amount of 
energy that is carried away by the spectator nuclei of the molecular target. In any case the 
models proposed here turn out to provide a simple and effective way of describing fusion 
reactions between hydrogen-like atoms. 
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VI. CONCLUSIONS 

This work proposes a simple and efficient model for the study of the screening enhancing 
effect on low-energy nuclear fusion reactions. In that model, the fusing atoms are considered 
hydrogen-like atoms whose electron probability density is used in Poisson's equation in order 
to derive the corresponding screened Coulomb potential energy. This way atomic excitations 
and deformations of the reaction participants can be taken into account. The derived mean- 
field potentials are then treated semiclassically, by means of the WKB, in order to derive 
the screening enhancement factor which is shown to be compatible with the experimentally 
obtained one for the H 2 (d, p) H 3 reaction, although some ambiguity remains regarding the 
molecular nature of the deuteron target. Moreover, by means of the proposed model the 
effect of a superstrong magnetic field on laboratory Hydrogen fusion reactions is investigated 
for the first time showing that despite the remarkable increase in the cross section of the dd 
reaction, the pp reaction is still too slow to justify experimentation. 
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FIGURE CAPTIONS 

Figure 1. The screening (acceleration) factor /i s with respect to the relative interaction 
energy of two fusing protons for various superstrong magnetic fields (in units of 10 12 G) 

Figure 2. The screening (acceleration) factor / ls with respect to the relative interaction 
energy of two fusing deuterons for various superstrong magnetic fields (in units of 10 12 G) 

Figure 3. The H 2 (d,p)H 3 astrophysical factor S (E) measured in keV-b with re- 
spect to the center of mass interaction energy E cm (keV) . The data (squares) are taken 



from Ref. ||20|| . The solid curve represents Eq. (|33|) , which makes use of the asymptotic 
form given in Ref. ||21|| . The dashed curve represents Eq. ( |3~4| ) where all four parameters 
S (0) , S (0) , S (0) , U e are fitted simultaneously. The dotted curve is obtained by adopting 
as a screening energy the value given by Eq. (|25| ) , while the dash-dotted curve stands for 
the astrophysical factor obtained by using Eq. (pTF) . 
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